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We study gravitational- wave emission from the quasi-circular, extreme mass ratio inspiral of com- 
pact objects of mass mo into massive objects of mass M ^ mo whose external metric is identical to 
the Schwarzschild metric, except for the absence of an event horizon. To be specific we consider one 
of the simplest realizations of such an object; a nonrotating thin-shell gravastar. The power radiated 
in gravitational waves during the inspiral shows distinctive peaks corresponding to the excitation 
of the polar oscillation modes of the gravastar. For ultra-compact gravastars the frequency of these 

■ peaks depends mildly on the gravastar compactness. For masses M ~ IO^Mq the peaks typically 
' lie within the optimal sensitivity bandwidth of LISA, potentially providing a unique signature of 
, the horizonless nature of the central object. For relatively modest values of the gravastar compact- 

■ ness the radiated power has even more peculiar features, carrying the signature of the microscopic 
properties of the physical surface replacing the event horizon. 

oo ' 

CJ ' I. INTRODUCTION 

cr. 

, One of the most elusive properties characterizing the black hole (BH) solutions of general relativity is the presence 
of their event horizon. Traditional electromagnetic astronomy can at best yield lower limits on the gravitational 
redshift corresponding to hypothetical surfaces replacing the event horizon (see for different viewpoints on this 
delicate issue). However, present and planned gravitational- wave (GW) detectors offer new prospects for "directly" 
observing BHs and probing their structure Q. 
"J ■ From a gravitational point of view the structure of compact objects is encoded in their free oscillation spectrum, 

■ i.e. in their quasinormal modes (QNMs) [EH- In an effort to point out the peculiar features of compact objects 
, whose external metric is identical to the Schwarzschild metric, but which do not possess an event horizon, in Ref. Q 

(henceforth Paper I) we studied the free oscillations of one of the simplest ultra-compact horizonless objects: a 
nonrotating thin-shell gravastar. Our analysis completed previous investigations by considering the thin shell 

as a dynamical entity. The QNM spectrum of a thin-shell gravastar is complex and profoundly different from that of 
a BH, mainly because of the different boundary conditions at the surface replacing the event horizon. 

As first proposed by Ryan [ll|, [l2|, an exquisite map of BH spacetimes can be constructed by observing the 
gravitational waveform emitted when a small compact object spirals into the putative massive BH at the center of a 
galaxy with the Laser Interferometer Space Antenna (LISA). Li and Lovelace refined the analysis showing that the 
' small object's tidal coupling encodes additional information about the metric of the spacetime [l^. In this paper 
we use the formalism developed in Paper I to show that any surface replacing the BH event horizon will produce a 
very specific signature in the gravitational signal emitted by the orbiting object because of the resonant scattering of 
gravitational radiation, which can be traced back to the different QNM spectrum of the two objects. In fact, here we 
show that the QNMs of ultra-compact thin-shell gravastars can be excited during the inspiral, whereas Schwarzschild 
QNMs can only be excited by particles plunging into the BH (see Refs. |14| - [l7| for a discussion of the analogous 
problem of particles orbiting neutron stars) . 

This work is very similar in spirit to a previous study by Kcsden and collaborators [T8| . There are two main 
differences between our work and theirs. The first difference is that Kesden et al. considered boson stars rather than 
gravastars as BH strawmen, so no "hard surface" replaces the event horizon in their case. The second difference is that 
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we compute the radiation in a consistent perturbative framework, instead of using "kludge" waveforms that become 
increasingly inaccurate in the relativistic regime. In this sense, this paper is the first "strong-field" calculation of the 
potential gravitational signatures of inspirals into horizonless objects. 

The plan of the paper is as follows. In Sec. Ill Al we review the equations describing axial and polar perturbations in 
the interior of a gravastar and the matching conditions with the ordinary perturbations of a Schwarzschild metric in 
the exterior spacetime. In Sec. IIIBl we summarize the perturbed Einstein equations in the Bardeen-Press-Teukolsky 
(BPT) formalism and we write down the source term appearing on the right-hand side of the BPT equation for 
orbiting pointlikc particles of mass mo ^ M. Then we discuss how the perturbation functions outside the shell (as 
obtained by the matching conditions derived in Paper I) can be used to solve numerically the BPT equation with a 
source given by the orbiting particle and to compute the radiated power. In Sec. lIIII we compare numerical calculations 
of the power radiated by BHs and different gravastar models and we stress potentially observable GW signatures of 
horizonless ultra-compact objects. We conclude by discussing possible extensions of our work. 



II. GRAVITATIONAL PERTURBATIONS BY A POINT PARTICLE 



The emission of gravitational waves by an extreme mass ratio binary system can be computed by a perturbative 
approach. One of the two compact objects is assumed to be an extended body of mass M , whose equilibrium structure 
is described by an exact solution of the Einstein equations. The second object is regarded as a point particle of mass 
jTiQ ^ M perturbing the gravitational field of the companion. This method has been applied to study gravitational 
radiation from particles orbiting BHs (see e.g. [l9l-l2ll|) and neutron stars [l3 - [l6l |. Here we apply the same perturbative 
approach to compute the power radiated by particles orbiting a thin-shell gravastar. 



A. Gravitational perturbations of the internal structure 



In this and the following sections we shall use the same notation as in Paper I. We consider a static thin-shell 
gravastar with metric [13, [IS] 
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where M is the gravastar mass and p = 3M/(47ra^) is the "energy density" of the interior region. The junction 
conditions require the induced metric to be continuous across the shell at r = a. This implies that /(r) must also be 
continuous at r ~ a, i.e. a = 1. In order to compute the radiation emitted by a gravastar perturbed by a massive 
point particle we must compute the gravitational perturbations both inside and outside the gravastar. Perturbations 
in the interior have been discussed in Paper I using the Regge- Wheeler (RW) gauge. Here we only recall some results 
that are useful to compute the power radiated by orbiting particles. 

In the de Sitter interior, both axial and polar perturbations can be reduced to the study of the master equation 
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where the tortoise coordinate is defined by dr/di\ ~ f{r) and /(r) is given by Eq. 
solution at the center (r = 0) is 



(2.3) 

with r < a. The regular 
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where C = {2M/a)^ = and F{a,b,c, z) is the hypergeometric function [23|. From Vt'" and its derivative we can 
obtain the Zerilli and RW perturbation functions as explained in Paper I. 
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For a thin-shell gravastar the background surface energy vanishes, E = 0, but the surface stress-energy tensor 
Q is, in general, nonvanishing. This implies that the perturbation functions are discontinuous across the shell. In 
Paper I we derived the matching conditions relating interior and exterior perturbations in the RW gauge. For axial 
perturbations these matching conditions read 



[N] = o, 



\/hhi 



0, 



(2.5) 



where [[. . .]] denotes the "jump" of a given quantity across the shell, i.e. the difference between the limits of the 
corresponding quantity as r — > a±. For a thin-shell gravastar Eqs. (j2.5p imply continuity of the RW function and its 
derivative across the shell. The treatment of polar perturbations is more involved and it yields the following relations 
for the jump of the polar metric functions across the shell: 
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The parameter Vs depends on the equation of state (EOS) on the thin shell, 9 = 6(S): 

9S 



(2.7) 
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and it has the dimensions of a velocity. A microscopic model of matter on the thin shell is needed for a microphysical 
interpretation of Vg , but (roughly speaking) this parameter is related to the sound speed on the shell. Wc shall follow 
the treatment in Paper I and treat Vs as a free parameter, although we will primarily focus on the (presumably more 
physical) range < < 1. 



B. The source term and the BPT formalism 



A detailed treatment of the perturbative approach used to compute the gravitational emission by a particle orbiting 
a polytropic neutron star can be found in Refs. [TsI - ITtI . [23|. Here we review the method with an emphasis on the 
modifications required to deal with thin-shell gravastars. 

The radial part >I';m(w,r) of the perturbation of the Weyl scalar S'i>4 is defined as 

^Uu;, r)^^Jdndt -2SU(^, 4>) [r^ 5^A{t, r, 0, 0)] e'^* , (2.8) 

where -2Sim{d,4') is a spin-weighted spherical harmonic of spin —2. The function ()2.8p can be expressed in terms of 
the Zerilli and RW perturbation functions {Zi{uj,r) and Yi{uj,r), respectively) as follows: 

^Uu^,r) = !:^^^^^±^[v--Yi + {W-- + 2iu;)A+Yi] (2.9) 



rV"(" + l) ^ypoi^^ ^ ^^poi ^ 2z^) A+Zi] 



4 

where 2n = (Z — + 2), A_|_ = d/dr^ + iu; = r^^Ad/dr + iuj and 

W^"" = ^(r-3M), (2.10) 

(nr + SAl ) 

The functions yP°' and are the well-known Zerilli and RW potentials: 
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with A = ?(/ + l)/2-l. 

The radial part of 6'^ 4 outside the shell can be computed from the matching conditions discussed in Section III Al 
and used as a boundary condition for the integration of the inhomogeneous BPT equation [lH, [26j 
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' — 2AIr and the source term Ti,„(a;, r) describes the point mass mg moving on a given orbit around the 
gravastar. In Ref. [l^ the solution of Eq. p.l4p is constructed in the general case of elliptic orbits. Eccentricity is 
expected to play an important role in extreme mass ratio inspirals [27I [2^. However, in the remainder of this paper 
we focus on circular inspirals. This simplifies our study and it is sufficient to prove our main point: the gravitational 
radiation from extreme mass ratio inspirals around horizonless objects is drastically different from the BH case. We 
mention in passing that our numerical code is capable of handling eccentric orbits, and the extension of our study to 
eccentric inspirals could be an interesting topic for future research. 

We further simplify the problem by using the so-called adiabatic approximation (i.e. we assume that the radiation 
reaction timescale is much longer than the orbital timescale). Under this assumption the trajectory of the particle is 
described by the geodesic equations for a mass mo moving on a circular orbit of radius Rq: 
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where t he dot indicates differentiation with respect to proper time. E is the energy per unit mass of the particle and 
wk = \/ M/Rq denotes the Keplerian orbital frequency. The source term can be written as 

,7r 
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where the functions sUim are explicitly given in Refs. [l6l[23|. 

The solution of Eq. (j2.14|) satisfying the boundary conditions of pure outgoing radiation at radial infinity and 
matching continuously with the interior solution can be found by the Green's functions technique. The amplitude of 
the wave at radial infinity can be shown to be [l5| 
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where Wim{<-^) is the Wronskian of the two independent solutions of the homogeneous BPT equation 
The two solutions ^I*;^ and ^^^^^ satisfy different boundary conditions: 
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Here Cbpt is the differential operator on the left-hand side of the BPT equation (|2.14p and ^'(^(Wja) is the radial 
perturbation of the Weyl scalar, constructed according to Eq. (|2.9p in terms of the perturbed metric functions in the 
interior and evaluated at the (exterior) surface of the gravastar. The integral in Eq. (j2.17p can be written in terms of 
and its derivatives [l^. In Eq. (|2.17p it is convenient to isolate the contribution of the Dirac 6 function: 
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Then the time-averaged energy-flux 
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can be written in terms of Ai„i{uj) as follows: 
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In order to evaluate ^'^^^ and we integrate the BPT equation by an adaptive Runge-Kutta method. Close to a 
resonance the solutions must be computed very accurately, since the Wronskian (|2.18p is the difference between two 
terms that almost cancel each other. When required, the tolerance parameter in the adaptive integration routines is 
decreased to achieve convergence. Since the orbital frequency is related to the orbital velocity v and to the semilatus 
rectum (which for circular orbits is simply p = Rq/M) by the relations 

V ^ (Mi^Kf' ^ p-'^\ (2.23) 

the energy flux can also be considered as a function of v or p. In the following we shall normalize E^ to the 
Newtonian quadrupole energy flux 

6" ^ (2,24) 

Then the energy flux emitted in gravitational waves normalized to the Newtonian quadrupole energy flux is given by 

Im 

The normalized energy flux p.25p can be computed up to w < 1/V6 ~ 0.408, which corresponds to the innermost 
stable circular orbit (ISCO) at Rq = 6M. The post-Newtonian expansion of the e nerg y flux P{v) for particles in 
circular orbit around Schwarzschild BHs has been studied by several authors [ll], [2^ |30|. The instability of circular 
orbits with Rq < 6M sets an upper bound on the velocity of the point mass. If the radius of the gravastar is larger 
than the ISCO (this typically occurs for fj, < 0.1666) the upper limit in v will be smaller. 

From the analytical form of the stress-energy tensor (|2.16p it is easy to see that, for each assigned /, a mode of the 
star is excited when the orbital frequency satisfies the resonant condition 

mujK = WQNM , (2.26) 

where wqnm is the QNM frequency. Thus we expect sharp peaks to appear at the values of v corresponding to the 
excitation of the gravastar QNMs for different values of the angular momentum parameter I. This offers an intriguing 
signature of the absence of event horizons, since the emitted power for a Schwarzschild BH does not show any peak. 
In fact one can easily check that the frequency of the fundamental QNM of a Schwarzschild BH is higher than the 
critical value ttjwk corresponding to a particle at the ISCO [IJ]. In other words, Schwarzschild QNMs can only be 
excited by particles plunging into the BH, while the QNMs of a gravastar can be excited during the inspiral. In the 
following Section we will compare the power emitted by a circular inspiral around a thin-shell gravastar to the power 
emitted by a circular inspiral around a Schwarzschild BH. 



III. GRAVITATIONAL FLUX FROM GRAVASTARS AND BLACK HOLES 



Thin-shell gravastar models are specified by two parameters: the gravastar compactness fj, = M/a and the sound 
speed parameter Vg that characterizes the EOS on the shell. Thin-shell gravastars are only one of the several possible 
models that can be explored (see e.g. [l^, [13) but we expect the qualitative results of our analysis to apply quite 
in general. The reason is that the main difference between gravastars and BHs comes from the different boundary 
conditions at the "surface" replacing the BH event horizon, rather than from the specific nature of this surface. 
Furthermore, as discussed below, peaks in the energy flux are more sensitive to the "global" properties of the gravastar 
(as determined by the compactness parameter /x) than to the microphysical model determining the matter distribution 
on the shell (which in our simplified case reduces to the specification of a value for Vs). Our numerical study covers 
the whole range in compactness (0 < /i < 0.5). We mainly focused on the most physical range of the EOS parameter 
(0 < 17^ < 1) but we also studied the supcrluminal case (w^ > 1), and we even allowed for models with < [3l|. 

The gravitational emission of a Schwarzschild BH perturbed b y a particle has been studied analytically and numeri- 
cally in great detail for both circular and eccentric orbits [T9l - [2ll [3^. Our purpose here is to compare and contrast the 
energy flux from particles orbiting Schwarzschild BHs to the energy flux from particles orbiting thin-shell gravastars. 
For each value of the gravastar parameters (/.t, w^) we integrate the perturbations equations (as described in Section HI]) 
for a point-like object of mass tuq moving on a circular orbit of radius Rq with orbital velocity v and we compute the 
energy flux (|2.25p . Our numerical work uses a modified version of the BPT code described in Ref. [3. The results 
obtained by the BPT formalism were verified using an independent code that integrates the Zerilli and Reggc- Wheeler 
equations. A slight variant of these codes was used to compute the flux from a particle orbiting Schwarzschild BHs. 
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The results are consistent with Refs. [l^ [2l| within an accuracy of about one part in 10^ (see Ref. [s^l for more 
details). 

From the results of Paper I, in the Schwarzschild limit /i 0.5 the real part of the QNM frequency tends to zero 
and to a very good approximation it is independent of Vs- For example, for ^ = 0.49999 and I = 2 we varied 
in the range [—2,2] in steps of Svg = 0.1 and we found that the real part of the modes is a constant within a part 
in 10^ {ujR = 0.235932), while the imaginary part has tiny variations in the range between w/ = 4.20 x 10^^ and 
U!i = 4.17 X 10~^. In order for a QNM to be excited by particles in circular orbits, the QNM frequency must be small 
enough to allow for the resonant condition (|2.26p . 
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FIG. 1: Dominant {I = 2) contribution to the energy flux for very high compactness and Vg 
)j, ~ 0.5 resonances are almost independent on v^). 
0.49997 , 0.49998 , 0.49999 , 0.499995 , 0.499999, respectively. 



0.1 (but when 

From right to left the resonant peaks correspond to /x = 



Figure [T] shows the dominant (/ = 2) contribution to the energy flux for gravastars with very high compactness. 
The frequencies of the lowest QNMs of a Schwarzschild BH are higher than those of an ultra-compact gravastar, 
and cannot be excited by particles in stable circular orbits. For this reason the power emitted by a BH (on the 
scale of this plot) is roughly constant. Resonance peaks do appear for gravastars, as expected, when wqnm = 2wk. 
Notice that these resonances are extremely narrow and they would get even narrower for I > 2. This is because the 
imaginary part of the excited modes is extremely small {2Mlui ~ 10""^, 10~^^ for I = 2 and I = 3 respectively) in the 
high-compactness limit /i — )- 0.5, leading to a corresponding decrease in the quality factor of the oscillations. Whether 
these resonances are actually detectable is an interesting question for LISA data analysis. The answer depends on 
dissipative mechanisms (besides gravitational radiation damping) that could affect the timescale of the oscillations, 
especially in the non-linear regime: see e.g. [itI |33| for discussions of this problem in the context of neutron star 
binary detection by Earth-based GW interferometers. 

Quite interestingly, gravastars that "try harder" to look like a BH (in the sense that their shell is closer to the 
Schwarzschild event horizon) are those for which the peak in the energy flux appears for smaller values of Table U 
lists the expected excited modes for different values of ^ corresponding to ultra-compact gravastars. 

One may worry that the resonance will eventually get out of the LISA band for gravastars having fi extremely close 
to the Schwarzschild value. The following naive argument suggests that this is not the case. The "thick shell gravastar" 
model by Mazur and Mottola predicts a microscopic but finite shell thickness £ ~ a/Ipits — 3 x 1O~^*(M/M0)^''^ cm, 
where Lpi is the Planck scale and rs is the Schwarzschild radius, so that the energy density and pressure in the shell 
are far below Planckian and the geometry can still be described reliably by Einstein's equations [13] . Our simplified 
model does not allow for a finite thickness of the shell, and a microscopic model of finite shells is required for a 
careful analysis of this problem. However, for the sake of argument, let us consider e = l/2 — /i— >0asa "thickness 
parameter" describing how far the gravastar shell can be relative to the BH horizon. A power-law fit of the QNMs of 
a thin-shell gravastar in the limit e — )■ yields /gw ~ 3.828 x (e x io~5")0.io73^ r^^^^ lower frequency sensitivity limit 
for LISA is dictated by acceleration noise. Assuming lower frequency cutoffs of /low = 10"^, 3 x 10^'^, 10~^, we find 
that the peaks will sweep out of the LISA band when e = 9.6 x 10~'^^,2.7 x 10~^^,2.0 x 10"^^, respectively. This 
estimate of the "minimum measurable deviation from a BH" is admittedly very sensitive to the fitting function we 
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TABLE I: Values of the compactness /i, angular momentum number /, QNM frequency, orbital velocity v and GW frequency 
vgw of the circular orbits which would excite the fundamental QNM of the gravastar for the given multipole. The Keplerian 
frequency is given in mHz and rescaled to a gravastar mass Afg ~ W^Mq. 





I McjQNM 
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{Me,/M)uGw (mHz) 


0.49997 
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0.1339 


0.4061 


4.328 
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0.1508 


0.3691 


4.873 


0.49998 
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0.1276 


0.3996 
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0.1429 


0.3625 


4.616 


0.49999 
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0.1180 


0.3893 
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0.1310 


0.3521 


4.232 


0.499995 


2 


0.1096 


0.3799 


3.543 


0.499999 


2 


0.0941 


0.3610 


3.041 



use and it may change when one considers thick shell gravastars, but it suggests that LISA has the potential to reveal 
solid surfaces replacing horizons even when these solid surfaces are very close to the location of the Schwarzschild 
horizon. 




FIG. 2: Left: The energy flux (summed up to / = 6) of GWs emitted by a small mass orbiting thin-shell gravastars with 
v1 = 0.1 and different values of /x (plotted as a function of the particle orbital velocity v) is compared with the flux for a 
Schwarzschild BH. All peaks (with the exception of the last two peaks on the right) are due to the excitation of QNMs with 
I = m. Right: same for vl = 0.1 and selected values of G [0.29, 0.49]. No QNMs are excited in this range. 



A relevant question is whether massive horizonless objects which are compact by the standard of (say) main sequence 
stars, but "only" as compact as neutron stars, can leave a signature on the gravitational signal emitted by small, 
inspiralling compact objects. In Fig. [2] we plot the normalized energy flux P{v) as a function of the orbital velocity 
for gravastar models with vl = 0.1 and compactness in the range 0.1 ^ ^ 0.49, as well as for a Schwarzschild 
BH. The total flux was computed by adding all multipoles (|m| < I) and by truncating the multipolar expansion at 
Z = 6. As discussed in Refs. [l^ HI, HO, [13 > ^ multipole of order I contributes to the total power as a correction of 
order p^~^ . Roughly speaking, a truncation at Z = 6 produces a relative error (in the non- resonant regime) of order 
= v^^ (but see [30| for a more careful discussion of the convergence properties of the post- Newtonian series). 
When fjL > 0.166 the ISCO is located outside the gravastar and we plot the energy fltix up to the ISCO velocity 
uisco — 0.408 (corresponding to Rq = 6M). For less compact gravastars, plots of the energy flux are truncated at 
the velocity corresponding to the location of the shell. 
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FIG. 3: Real part (left) and imaginary part (right) of gravastar QNMs with / = 2 as a function of compactness for several fixed 
values of (as indicated in the legend). For clarity in illustrating the "selection rules" that determine QNM excitation during 
inspiral we only show the weakly damped part of the QNM spectrum (compare Fig. 2 and Fig. 6 in Paper I). For > 0.8 
the real part of the frequency is plotted down to the critical minimum compactness at which the imaginary part crosses zero 
within our numerical accuracy, the horizontal line at 2Mtj_R ~ 0.2722 corresponds to twice the orbital frequency of a particle 
in circular orbit at the ISCO: only QNMs below this line can be excited during a quasi-circular inspiral. 




J ■ ■ I i_ 



FIG. 4: Real part (left) and imaginary part (right) of gravastar QNMs with Z = 3 as a function of compactness for several fixed 
values of (as indicated in the legend). For clarity in illustrating the "selection rules" that determine QNM excitation during 
inspiral we only show the weakly damped part of the QNM spectrum. In the left panel, the horizontal lines at 2Mu}r ~ 0.1361 
(2Mtj_B ~ 0.4082) correspond to the orbital frequency (or three times the orbital frequency) of a particle in circular orbit at 
the ISCO. Perturbations with / = 3,m = 1 can excite the QNMs below the first line, while perturbations with I = m = 2) can 
excite QNMs below the second line. 



The complex structure of the spectrum for values of ^ smaller than about 0.2 is best understood by considering 
the real and imaginary parts of the weakly damped QNM frequencies of a gravastar (see Fig. [3]). For clarity in Fig. [3] 
we only plot weakly damped QNMs, but our general arguments apply also to the second, "ordinary" family of QNMs 
(cf. Paper I). In particular, from Fig. 2 and Fig. 6 in Paper I it should be clear that QNMs will be excited for low 
values of the compactness and when ^ is very close to the BH value ^ = 1/2. Besides these "ultracompact" modes, 
only QNMs whose real part lies helow the horizontal line in the left panel (corresponding to twice the ISCO orbital 
frequency for a particle in circular orbit) can be excited. 

Fig. ini clarifies that the range of ^ over which QNMs can be excited depends on Vg. For v1 ~ 0.1 (the case 



9 



considered to produce the energy fluxes of Fig. [2]) QNM frequencies that can be excited by resonant inspirals only 
exist for fj, < 0.21 (left panel of Fig. ^ or for /i > 0.49997, i.e. when the thin shell is extremely close to the location 
of the BH horizon (Fig. [T]). The real part of the corresponding QNM frequency has a local maximum at /i w 0.15. 
Correspondingly, the I = 2 QNM peak visible in the energy flux of Fig. [2] occurs later in the inspiral for the yu = 0.15 
model than it does for the /.i = 0.10 and yu = 0.20 models. 

In Fig. [2]the I = 2 and I = 3 peaks for yu = 0.20 are well separated in frequency and an "antiresonance" is visible to 
the right of the I = 2 resonance. The nature of this antiresonance can be explained by a simple harmonic oscillator 
model [l^. In the inset of the left panel of Fig. [2] we plot both the resonance and antiresonance as functions of 
the Keplerian orbital frequency of the particle AIluk for /.i = 0.2 and I = 2 (dashed green line). A fit using the 
simple harmonic oscillator model of Ref. [l^ (red line) reproduces the qualitative features of both resonance and 
antiresonance: in this specific case the fit gives 2AIujji ~ 0.07257 and 2Mu!j ~ 2 x 10^^, while QNM calculations 
using the resonance method yield 2Mujr ^ 0.07257 and 2Mui 4 x lO"*'. 

Modes with I > 2 are typically harder to excite because of their higher frequencies and lower quality factors. 
However, because of the complex "selection rules" illustrated in Fig.[3]for / = 2, sometimes only resonances with I > 2 
will be visible. When = 0.1 and fi > 0.21 only modes with I > 2 can be excited, and only narrow I = 3 resonances 
can be seen in Fig. [5] when the compactness /i = 0.25 (cf. Fig. d]). 

When I = 2 the imaginary part of one QNM with ~ 0.1 crosses zero within our numerical accuracy at the 
"critical" compactness /i ~ 0.21, possibly signaling a (marginal) nonradial instability of the gravastar, and no QNMs 
can be excited for 0.21 < n < 0.49997. In this compactness range the energy flux emitted by either the gravastar or 
the BH is mostly due to the orbital motion and it only depends on the compactness of the central object. The right 
panel of Fig. [5] shows that the flux emitted by a gravastar approaches the BH flux "from below" as the compactness 
increases. For fj, ~ 0.35 the gravastar flux is almost indistinguishable from the BH flux and for /i > 0.35 a gravastar 
radiates slightly more than a BH. This is due to the fact that the emitted power "feels" the contribution of resonances, 
which in this case correspond to orbits smaller than the ISCO but do nevertheless contribute to increase the slope 
of the curve. A similar trend can be seen in neutron star calculations in regions of the parameter space where the 
contribution from resonances is negligible [iGj . 

If gravastars or other horizonless objects have astrophysical reality, the presence or absence of resonant peaks in the 
GW flux can provide interesting information on the microscopic properties of the physical surface replacing the event 
horizon. Suppose for example that we can estimate the compactness of a massive object by independent means (e.g. 
by electromagnetic observations). Even within our simple thin-shell model, the range in frequency where resonances 
in the GW emission from EMRIs are allowed changes with v^. For example, ii = 0.1 resonances can exist when 
/i < 0.21 OT ^ 0.5, but if = 0.3 they can exist when ^ < 0.1, 0.33 < /x < 0.36 and ^ 0.5. Similar results also 
hold when ^ = 3, as shown in Fig. For example, if v"^ = 0.1 resonances can exist when /i < 0.27 or jj. ~ 0.5, but 
if = 0.3 they can exist when < 0.2, 0.37 < ^ < 0.41 and ^ ~ 0.5. So, in general, the range of fi where QNM 
frequencies can be excited by a circular inspiral depend on the value of I. In the Schwarzschild limit {fj, ~ 0.5) QNM 
frequencies are excited for any I, but higher-^ modes have a tiny imaginary part {2AIujj ^ lO^^'' for I = 3) and they 
are more difficult to detect than the dominant {I = 2) modes. 

If an EMRI is detected, the existence of these selection rules (in the form of compactness regions where resonances 
can or cannot exist) in principle allows for null tests of the existence of an event horizon for objects of the given 
compactness. Similar arguments can presumably be made for more complex (or contrived) gravastar models. 



IV. CONCLUSIONS 



We have shown that gravitational radiation from EMRIs can be used to tell the presence or absence of an event 
horizon in a compact, massive object. More specifically, we have shown that the resonant excitation of the oscillation 
modes of a gravastar in the LISA band is a potentially observable signature of the surface replacing the event horizon. 
For thin-shell gravastar models there is a range of compactness (e.g., fi < 0.21 and /i > 0.49997 for = 0.1) where 
this resonant scattering can occur. 

More detailed data analysis studies (possibly including refined microphysical models of this "solid surface") are 
necessary to determine the detectability of resonant peaks, especially for ultra-compact gravastars. The extension of 
our results to rotating gravastar models presents a challenge because of the difficulties in finding plausible rotating 
gravastar solutions and because of the ergoregion instability that affects some rotating gravastar models [ssl - fSTt . 
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